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ABSTRACT 

Longwave  Infrared  (LWIR)  radiation  eomprising  atmospherie  and  surfaee  emissions 
provides  information  for  a  number  of  applications  including  atmospheric  profiling,  surface 
temperature  and  emissivity  estimation,  and  eloud  depietion  and  eharaeterization.  The  LWIR 
spectrum  also  contains  absorption  lines  for  numerous  molecular  species  which  can  be  utilized 
in  quantifying  speeies  amounts.  Modeling  the  absorption  and  emission  from  gaseous  speeies 
using  various  radiative  transfer  codes  such  as  MODTRAN-4  and  FASE  (a  follow-on  to  the 
line-by-line  radiative  transfer  eode  FASCODE^)  provides  insight  into  the  radiative  signature 
of  these  elements  as  viewed  from  an  airborne  or  spaee-borne  platform  and  provides  a  basis 
for  analysis  of  LWIR  hyperspeetral  measurements. 

In  this  study,  a  model  platform  was  developed  for  the  investigation  of  the  passive  outgoing 
radiance  from  a  scene  eontaining  an  effluent  plume  layer.  The  effeets  of  various  scene  and 
model  parameters  including  ambient  and  plume  temperatures,  plume  eoneentration,  as  well 
as  the  surfaee  temperature  and  emissivity  on  the  outgoing  radiance  were  estimated.  A  simple 
formula  relating  the  various  components  of  the  outgoing  radiance  was  used  to  study  the  scale 
of  the  component  contributions.  A  number  of  examples  were  given  depicting  the  spectral 
radiance  from  plumes  eomposed  of  single  or  multiple  effluent  gases  as  would  be  observed  by 
typieal  airborne  sensors.  The  issue  of  deteetability  and  speetral  identifieation  was  also 
diseussed. 

1.  INTRODUCTION 

Detection  of  effluents  from  remotely  based  platforms  provides  a  meehanism  to  monitor  a 
variety  of  environmental  and  geologic  conditions  over  a  large  region.  Volcanic  eruptions 
produce  large  amounts  of  SO2,  H2O  and  CO2  whieh  ean  be  detected  from  high  altitude 
sensors.  Oil  spills  can  be  monitored  by  their  speetral  signature.  There  is  a  strong  interest  in 
the  military  to  improve  our  ability  to  deteet  and  identify  biological  and  chemical  releases  in 
the  battlespace  arena.  For  example,  the  Multispectral  Thermal  Imager  (MTI)  sensor  was 
designed  to  aid  in  the  deteetion  of  the  proliferation  of  nuclear  weapons"^. 


*  This  work  was  sponsored  by  the  Department  of  Air  Force  under  Contract  F19628-00-C-0002.  Opinions, 
interpretations,  conclusions  and  recommendations  are  those  of  the  authors  and  not  necessarily  endorsed  by  the 
United  States  Air  Force. 
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The  ability  to  model  hyperspeetral  seenes  whieh  eontain  effluent  emissions  from  a  staek  or 
other  souree  provides  the  initial  foundation  for  further  efforts  in  effluent  deteetion  and 
remote  identification.  The  simulation  of  the  scene  radiance  for  plume  cases  was  primarily 
used  for  development  of  applications  for  a  hyperspeetral  system  analysis  model  termed  the 
Forecasting  and  Analysis  of  Spectroradiometric  System  Performance  (FASSP^).  It  provides  a 
means  for  understanding  the  effects  of  plumes  on  the  outgoing  radiance  field  and  provides  a 
testbed  for  the  ongoing  effluent  detection  and  identification  study. 

In  this  study,  FASE  is  used  to  produce  estimates  of  the  atmospheric  spectral  opacity  (layered 
optical  depth)  from  a  scene  containing  an  effluent  plume  layer,  which  when  used  in 
conjunction  with  a  simple  radiative  transfer  algorithm  can  produce  an  estimate  of  the  at- 
sensor  radiance.  In  a  previous  study,  the  effects  of  various  scene  and  model  parameters 
including  surface  and  ambient  temperature,  plume  temperature  and  concentration  and 
spectral  bandwidth  on  the  outgoing  radiance  were  estimated^.  Results  of  trade  studies  on  the 
magnitude  of  the  various  parameter  contributions  to  the  outgoing  radiance  showed  that  plume 
concentration  and  temperature  were  the  most  significant  contributors.  In  this  paper,  the 
relationship  between  these  parameters  will  be  explored  in  more  detail.  Using  gaseous 
absorbances  from  a  commercially  available  spectral  library,  at-sensor  radiances  can  be 
modeled  for  combinations  of  gases  not  found  in  standard  line-by-line  spectral  databases  such 
as  HITRAN-96^.  This  provides  the  opportunity  to  simulate  the  measured  signal  from  a 
controlled  release  of  gases  as  observed  from  an  airborne  or  spacebome  sensor.  Some 
examples  of  these  simulations  are  given  in  later  sections  of  this  paper.  Using  the  FASE  line- 
by-line  radiative  transfer  code,  estimates  of  the  at-sensor  radiance  can  be  made  for  a  selection 
of  sensor  bandwidths,  effluent  concentrations  and  temperatures  and  scene  backgrounds.  An 
example  of  the  signal  from  an  atmosphere  with  a  layer  of  enhanced  SO2  is  shown  in  Eig.  la. 
The  difficulty  in  identifying  the  effluent  solely  on  its  radiance  is  clearly  seen.  However,  if  an 
estimate  of  the  background  radiance  can  be  obtained,  then  the  radiance  difference  due  to  the 
enhanced  effluent  can  be  estimated.  The  characteristic  double  peaked  SO2  absorption  signal 
in  the  8.5  |im  region  is  easily  recognized  in  Figure  lb  where  the  spectral  range  of  the  plots 
has  been  narrowed. 

2.  RADIATIVE  TRANSFER  FOR  PLUME  SCENARIO 

To  model  the  at-sensor  radiance  for  a  scene  containing  a  plume,  some  modeling  assumptions 
were  made:  1)  the  enhanced  effluent  amounts  reside  wholly  within  a  thin  plume  layer 
(nominally  10  m),  2)  the  temperature  is  constant  over  the  entire  plume  layer,  and  3)  the 
surface  conditions  (temperature  and  emissivity)  for  the  background  and  underlying  the  plume 
are  identical.  Figure  2  shows  a  diagram  of  the  various  components  of  the  background  and 
plume  radiance  computation.  Ideally,  the  signature  of  the  effluent  can  be  obtained  by 
subtracting  the  radiance  of  a  nearby  effluent-free  location  from  the  radiance  through  the 
plume.  Radiative  transfer  equations  for  both  the  background  and  through-plume  cases  are 
presented  below. 


2,1  Clear-sky  Radiance 


The  clear  or  background  radiance  can  be  found  by  summing  the  contributions  from  three 
components:  1)  the  path  radiance  Lpath,  2)  the  surface  emitted  radiance  transmitted  to  the 
sensor  Ls/c,  and  3)  the  reflected  downwelling  radiance  transmitted  to  the  sensor  Lreji.  For  a 
single  wavelength  or  narrow-band  channel  the  equation  can  be  written  as, 

^hkg  ~  ^path  ^refl  ^sfc  ~  ^iip  ^ sfc  )^atm^dii  (  ^  ) 

where  s^c  and  Tg/c  are  the  emissivity  and  temperature  of  the  surface,  respectively,  B  is  the 
Planck  function,  Lup  and  Ldn  are  the  atmospheric  upwelling  and  downwelling  radiances. 
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Figure  1.  Plots  of  the  radiance  in  units  of  |lW  cm'  sr'  |im'  (jiflicks)  from  a  scene  with  and 
without  an  enhanced  SO2  concentration  layer.  The  upper  plots  provide  the  simulated 
radiances  at  1  cm''  resolution  while  the  lower  plots  represent  the  difference  in  the  two 
radiance  curves.  Curves  in  a)  represent  the  signal  across  the  broad  LWIR  spectral  band  while 
b)  narrows  the  spectral  range  to  that  which  covers  the  primary  SO2  absorption  band. 


respectively,  and  tam  is  the  atmospheric  transmittance  from  surface  to  sensor.  Each  of  the  full 
atmospheric  parameters  can  be  separated  into  individual  components  arising  from  below, 
above  and  within  the  plume  layer  as  depicted  in  Fig.  2.  For  the  background  case  there  is  no 
plume  so  the  plume  layer  is  treated  as  a  thin  model  layer  with  appropriate  amounts  of  all 
gases.  Therefore,  we  can  define  the  atmospheric  parameters  in  (1)  by  their  component  values, 

^ atm  ^ p*  ^ a 


where  the  sub/superscripts  b,  a  and  p*  on  the  transmission  and  radiance  represent  below, 
above  and  in-plume  layer  values,  respectively.  The  above-plume  transmission  ta  refers  to  the 
path  from  the  top  of  the  plume  layer  to  the  sensor.  This  also  applies  to  the  atmospheric 
transmission  tam  and  the  upwelling  radiance  Lup‘^.  The  downwelling  radiance  Ldn  is  measured 
over  the  entire  atmospheric  path  from  top-of-atmosphere  to  the  plume  layer.  Flpon 
substitution  (1)  becomes. 


{path) 
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Figure  2.  Graphic  depicting  the  contributions  from  the  surface,  plume  layer  and  atmosphere 
to  the  at-sensor  radiance.  The  model  levels  are  shown  as  a  function  of  altitude. 


It  can  be  shown  that  some  eomponents  provide  a  small  or  negligible  eontribution  to  the 
overall  signal.  The  plume  layer  and  below  plume  layer  transmissions  th  and  tp*_  respeetively, 
are  very  elose  to  1.  The  downwelling  radianee  from  below  the  plume  layer  dn  is  mueh 
smaller  than  that  from  above  the  plume  layer.  Likewise  for  the  upwelling  radianee  l!^up-  If 
these  assumptions  are  followed  then  (2)  reduees  to  Lup  =  L^up,  Ldn  =  L^^dn  and  tam  =  L  and  (1) 
ean  be  written  in  the  form, 

^bkg  ~  ^ path  ^refl  ^ sfc  ~  ^up  ~  ^ sfc  )  K^dn  ^ sfc^  sfc  \  K  '  (4) 

The  atmospherie  parameters  L^^up,  L^dn  and  ta,  ean  be  estimated  from  initial  seene  eonditions 
using  an  appropriate  radiative  transfer  model.  The  surfaee  temperature  Ts/c,  may  be  known 
from  supplementary  data,  while  the  surfaee  emissivity  £sfc,  a  funetion  of  wavelength,  must  be 
retrieved. 


2,2  Radiance  Through  a  Plume 


When  a  plume  layer  is  added  to  the  seene,  the  eontributions  to  eaeh  of  the  three  eomponents 
beeome  more  eomplex.  While  most  radiative  transfer  models  do  not  ealeulate  the  individual 
eomponents  assoeiated  with  the  inelusion  of  a  plume  (they  do  not  treat  the  plume  layer  any 
differently  from  one  with  standard  amounts  of  effluent),  it  is  of  interest  to  determine  the  seale 
of  the  eontributions  from  eaeh  of  the  loeations  denoted  by  a  solid  eirele  in  Fig.  2.  The 
radianee  equation  for  a  ease  with  a  plume  layer  has  the  form. 


r  =  If  If  If 

^ plume  ^ path  '  '  ^sfc 

=  V  +L^  t  t 

up  up  p  a 


{path) 

(refl) 

(sfc) 


(5) 


where  the  sub/superseripts  b  and  a  on  the  transmission  and  radianee  represent  below  and 
above-plume  layer  values,  respeetively,  and  where  the  sub/superseripts  p  on  the  transmission 
and  radianee  represent  the  in-plume  layer  value  while  p*  'm  (2)  above  represented  the  same 
layer  but  for  a  elear  seene.  In  this  ease  the  transmission  through  and  radianee  emitted  from 
the  plume  layer  ean  not  be  negleeted.  Some  terms  in  (5)  have  a  relatively  small  impaet  on  the 
overall  at-sensor  radianee.  Griffin,  et  al  ^  eomputed  the  pereent  eontribution  for  eaeh  of  the 
terms  in  (2)  and  (5)  for  a  seleetion  of  simulated  effluent  emissions.  The  dominant  term  for 
baekground  (no  plume)  eonditions  is  the  surfaee  emission  term  whieh  eontributes  to  93-95  % 
of  the  radianee.  Emission  from  the  plume  is  the  dominant  term  as  expeeted  in  the  effluent 
absorption  band.  The  results  showed  that  refleetion  terms  are  mostly  negligible  beeause  of 
the  low  refleetanee  of  the  surfaee  (generally  less  than  0.1  in  the  LWIR)  and  the  path  radianee 
(exeluding  the  plume  eontribution)  eontributes  5-8%  of  the  total  radianee. 


2,3  Differential  Radiance 


The  radiance  difference  AL  between  the  scene  with  a  plume  and  the  background  can  be  found 
by  subtracting  (2)  from  (5),  and  assuming  the  transmission  through  the  “clear”  thin  plume 
layer  tp*  is  1, 


(6) 


To  simplify  the  equation,  we  will  apply  the  assumptions  made  above  for  the  upwelling  and 
downwelling  radiance  below  the  plume  (L^up  «  L^up,  dn  «  dn)  and  for  the  transmittance 
for  that  layer  (4  ~  1)  to  (6),  which  now  reduces  to, 

Ai = <,  (i  -  b  k  1-  1+  (i  -  ){s  k  ] «,  -  (i + } } .  (7) 

The  third  component  within  the  parenthesis  is  a  form  of  the  reflected  downwelling  radiance. 
From  Griffin,  et  al  the  contribution  from  the  reflected  downwelling  radiance  accounts  for 
less  than  3  %  of  the  total  radiance  differential  value  depending  on  the  plume  temperature  and 
transmission  and  the  surface  emissivity.  Neglecting  the  downwelling  radiance  term  simplifies 
(7)  considerably  and  is  done  here.  The  equation  reduces  to  the  approximate  form, 

Ais(.(l-<,){skl-«-./.7SkJ}  .  (8) 

which  illustrates  the  three  factors  that  appear  to  have  the  most  effect  on  the  detection  of  an 
effluent:  1)  the  atmospheric  transmission  from  the  plume  to  the  sensor,  2)  the  emissivity  of 
the  plume,  and  3)  the  relative  radiance  difference  between  the  plume  and  the  surface.  The 
first  factor,  the  atmospheric  transmission  influences  the  effluent  signature  by  attenuating  the 
overall  signal.  It  affects  both  the  background  and  plume  emissions  equally  (assuming  the 
atmosphere  conditions  are  unchanging  over  the  scene).  The  second  factor,  the  emissivity  of 
the  plume  is  a  function  of  the  transmission  of  the  plume  layer  and  therefore  is  directly  related 
to  the  concentration  of  effluent  in  the  plume  layer.  Lower  plume  layer  transmission  (a  result 
of  increased  effluent  amount),  acts  to  increase  AL.  Conversely,  the  plume  emissivity  tends 
toward  zero  (as  does  AL)  when  the  effluent  concentration  is  small.  The  third  factor  which 
relates  the  plume  layer  radiance  and  the  surface  emitted  radiance  is  a  function  of  the  plume 
and  surface  temperature  and  the  surface  emissivity.  Here  the  contrast  in  the  plume  versus 
surface  temperature  affects  not  only  the  magnitude  of  the  signal  but  the  direction  of  the 
signal  as  well.  If  the  surface  temperature  is  significantly  greater  than  the  plume  temperature 
(i.e.,  hot  sand,  pavement  or  concrete),  AL  can  be  negative;  the  plume  would  look  cooler  than 
the  surrounding  objects  in  the  scene.  If  the  plume  and  surface  temperature  are  approximately 
equal  then  a  AL  near  zero  would  be  seen.  An  example  of  this  would  be  the  detection  of 
plumes  at  distances  away  from  the  source  where  the  plume  has  had  time  to  mix  with  the 


ambient  air.  Even  extremely  large  eoncentrations  of  effluent  might  be  undetectable  if  the 
temperature  contrast  were  not  present.  The  word  “might”  is  used  here  to  remind  us  that 
techniques  which  utilize  other  regions  of  the  spectral  band  (outside  of  the  LWIR)  may  be 
able  to  detect  either  an  emissive  temperature  difference  or  possibly  a  reflected  sunlight 
component  of  the  plume.  If  the  plume  temperature  is  greater  than  the  background 
temperature,  then  detection  of  the  effluent  is  possible.  In  fact,  both  the  effluent  concentration 
and  the  temperature  difference  work  to  enhance  the  signal:  as  concentrations  rise  the  plume 
transmission  decreases  and  the  second  factor  approaches  one. 

The  question  becomes  which  is  the  most  important  factor  in  modulating  the  cumulative 
radiance.  Changes  in  the  plume  temperature  are  directly  observed  in  the  Planck  emitted 
radiance,  and  similarly  for  the  surface  temperature.  Modifying  the  effluent  concentration 
results  in  an  attenuation  of  the  below  plume  radiance  most  significant  in  the  effluent 
absorption  band  region,  and  in  a  change  in  the  emissivity  characteristics  of  the  plume  (tp  ~  1- 
£p  in  the  LWIR).  The  increase  in  plume  emission  is  partially  balanced  by  a  decrease  in  the 
surface  emitted  radiance  passing  through  the  plume  and  reaching  the  sensor.  It  can  be  seen 
that  the  effluent  concentration  modulates  the  radiance  but  the  plume  temperature  determines 
both  the  baseline  radiance  and  whether  the  plume  scene  radiance  is  less  than  the  background 
radiance  (as  in  the  hot  surface  case)  or  greater  than  the  background  radiance  (warm  plume  or 
cool  surface).  In  the  case  where  the  plume,  ambient  and  surface  temperatures  are  identical,  it 
would  be  very  difficult  to  detect  an  effluent  from  the  radiance  signal  alone. 

An  estimate  of  the  effluent  concentration  can  be  obtained  from  the  plume  layer  transmission 
tp  by  solving  (10)  for  tp  , 


AL 

K  ~  ^sfc^  ^sfc  .  1 


(9) 


The  plume  layer  transmission  tp  =  10  '^  where  a  is  the  absorbance  (related  to  the  optical  depth 
by  a  constant  factor)  of  the  plume.  Therefore,  the  equation  for  a  has  the  form, 

a  =  logio  I X I  -  logio  I X  -  AL  I ,  where  x  =  t^{B  [l^  ]-  ] }  .  (10) 

In  (10)  it  is  clear  that  the  value  of  x  must  be  greater  than  AL.  As  AL  approaches  zero  (plume- 
free  case),  a  goes  to  zero  as  expected  and  when  zlL  =  x-1 ,  the  value  for  a  is  just  logio\x\.  It  is 
unlikely  the  latter  value  has  any  physical  significance;  its  occurrence  is  most  likely  a  factor 
of  the  assumptions  made  above.  As  AL  approaches  x,  a  becomes  increasingly  large.  To 
improve  the  accuracy  of  the  absorbance  estimate,  a  small  correction  factor  could  be  added  to 
X  to  include  the  effects  of  the  neglected  reflected  downwelling  component  in  (7). 

2,4  Detection  Issues 


The  detection  of  an  effluent  plume  from  an  airborne  or  spaceborne  platform  using  an  LWIR 
technique  similar  to  the  one  described  in  the  previous  sections,  requires  accurate  knowledge 


of  the  intervening  atmosphere  (spectral  transmission),  the  underlying  surface  (temperature 
and  emissivity),  and  the  spectral  absorption  characteristics  of  candidate  gaseous  species. 
Determination  of  the  background  radiation  is  the  most  important  piece  of  information  needed 
and  probably  the  most  difficult  to  estimate  accurately.  Errors  can  result  from  assuming  the 
radiance  from  a  nearby  clear  region  is  representative  of  the  radiance  underlying  the  plume.  If 
a  segment  of  the  spectral  band  is  known  to  be  free  of  effluent  absorption,  then  an  estimate  of 
the  effective  surface  radiance  Te  =  £sB[Ts]  might  be  derived.  However,  knowledge  of  the 
variation  in  emissivity  with  wavelength  is  needed  to  accurately  associate  the  derived  value  to 
other  spectral  bands.  The  separation  of  temperature  and  emissivity  is  a  fundamental  problem 
in  the  LWIR  and  several  techniques  have  been  published  to  address  this  issue^’^’^°’'\ 

The  above  approach  assumes  that  the  effluent  species  was  known  and  the  analyst  knew 
where  to  look  spectrally  to  find  the  telltale  signal.  Without  this  knowledge,  a  spectral 
matched  filter  technique  or  some  other  approach  must  be  used  to  provide  this  information. 
Using  a  mathematical  approach,  gaseous  absorption  signatures  from  a  spectral  library  are 
compared  with  the  measured  radiance  (actually  AL)  to  detect  and  identify  the  effluent 
species. 

3.  LWIR  RADIATIVE  TRANSFER  MODEL 

The  technique  used  in  this  study  to  model  the  spectral  radiance  from  a  scene  containing  a 
plume  layer  consists  of  the  following  steps:  1)  compute  the  layer  optical  depths  using  a  line- 
by-line  model  such  as  EASE  for  atmospheric  conditions  approximating  the  plume-free  scene 
conditions,  2)  merge  the  spectral  absorptivity  of  the  selected  effluents  into  the  model  plume 
layer,  and  3)  compute  the  atmospheric  at-sensor  radiance  components  using  (3)  and  (5).  The 
simulation  of  the  scene  radiance  for  plume  cases  was  primarily  used  for  the  development  of 
applications  for  FASSP^.  It  has  proven  a  useful  tool  for  the  understanding  of  effects  of 
plumes  on  the  outgoing  radiance  field  and  provides  a  testbed  for  the  ongoing  effluent 
detection  and  identification  study  .  A  brief  description  of  the  above  steps  follows. 

3.1  Plume  Scene  Modeling  Steps 

Using  atmospheric  properties  appropriate  to  the  desired  model  parameters,  EASE  was  run  to 
produce  optical  depths  at  user-defined  layers.  Since  only  atmospheric  parameters  were 
calculated,  no  input  surface  information  was  needed.  This  approach  allowed  the  definition  of 
surface  parameters  at  the  time  of  the  final  radiance  component  calculation,  without  the  need 
to  rerun  the  line-by-line  model.  The  full  resolution  layered  optical  depths  were  reduced  to  a 
resolution  of  0.1  cm'^  using  the  procedure  outlined  in  Griffin  et  al 

The  characterization  of  the  plume  layer  was  done  following  the  EASE  calculations  of  the 
layered  atmospheric  optical  depths.  The  effluent  dispersion  was  derived  from  a  gaussian 
plume  model.  This  provided  an  estimate  of  the  relative  amount  of  effluent  at  a  specified 
distance  away  from  the  source  (in  this  case  a  stack  of  known  height,  diameter  and  outflow 
parameters).  The  effluent  library  contained  the  absorptivity  as  a  function  of  wavelength  and 
integrated  gas  amount  (in  PPM-m)  for  each  effluent.  Using  the  known  effluent  emission  rate 


(in  g  s'^)  and  the  estimated  dispersed  gas  amount  from  the  plume  model,  the  eorresponding 
speetral  absorptivity  of  the  gas  to  plaee  into  the  plume  layer  was  determined.  In  faet,  multiple 
gases  eould  be  merged  together  using  individual  emission  rates  and  speetral  absorptivities  to 
obtain  a  plume  layer  representative  of  a  simultaneous  release  of  many  gases.  For  these  eases, 
it  was  assumed  that  a  single  temperature  was  representative  of  all  gases  in  the  plume  for  a 
speeifie  loeation  away  from  the  staek.  Also,  the  speetral  resolution  of  the  effluent 
absorbanees  differed  from  the  LWIR  model  speetral  resolution.  A  eubie  spline  interpolation 
seheme  was  employed  to  eonvert  the  absorbanees  to  a  standard  model  resolution  of  0.1  om'\ 

A  simple  radiative  transfer  (RT)  eode  was  developed  to  utilize  the  eomputed  layered 
transmittanee  information  along  with  known  atmospherie  temperature  profiles  and  modeled 
surfaee  temperature  and  emissivities  to  eompute  the  upwelling,  downwelling,  and  surfaee 
emission  eomponents  of  the  atmospherie  radiative  transfer.  The  assumption  of  a  non- 
seattering  (emission  and  absorption  only)  plane -parallel  atmosphere  was  used  to  simplify  the 
eomputations.  These  eomponents  ean  be  eomputed  for  both  the  with-plume  and  without- 
plume  atmospherie  eases. 

FASSP  allows  the  user  to  simulate  the  at-sensor  radianee  for  a  variety  of  hyperspeetral 
imagers  ineluding  VNIR/SWIR  sensors  sueh  as  AVIRIS^^  and  HYDICE^"^  as  well  as 
MWIR/LWIR  sensors  sueh  as  SEBASS'^.  As  a  final  step,  the  speeifie  sensor  response 
funetion  was  eonvolved  with  the  modeled  radianees  to  provide  an  estimate  of  the  sensor- 
speeifie  measured  radianee  for  the  plume  seenario. 

3.2  Modeling  Examples 

The  LWIR  RT  eode  was  used  in  eonjunetion  with  plume  dispersion  and  sensor  models  to 
estimate  the  at-sensor  radianee  from  a  seene  eontaining  a  gaseous  plume.  Figure  3a  depiets 
the  derived  spectral  radiances  (shown  as  equivalent  brightness  temperatures,  Tb)  for  a  scene 
containing  an  active  plume  composed  of  a  single  gas  (methanol)  emitted  from  a  stack  at  20 
m.  The  specific  parameters  for  this  scenario  are  given  in  Table  1.  Some  of  the  parameters  are 
used  in  the  plume  concentration  analysis  (stack  height,  diameter,  exit  velocity,  gas  emission 
rate,  wind  speed,  etc.)  while  others  are  used  in  the  radiative  transfer  equations  shown  above 
(sensor  altitude,  surface  and  gas  temperature)  and  in  the  sensor  definition  (sensor  type  and 
OPD).  Included  in  Fig.  3a  are  the  background  Tb  (dotted  line  in  Fig.  3a)  along  with  plots  of 
the  differential  Tb  (Fig.  3b),  the  plume  emissivity  (Fig.  3c)  and  the  effluent  absorbance  (Fig 
3d)  as  derived  from  (10).  While  the  calculation  of  the  absorbance  is  based  upon  radiance 
calculated  from  a  scene  containing  prescribed  gaseous  amounts,  it  does  still  provide  a 
consistency  check  on  the  methodology.  For  this  case,  the  plume  signature  is  easy  to  detect 
due  to  a  combination  of  enhanced  gas  concentration  and  a  warm  plume  temperature 
(approximately  20  K  higher  than  the  surrounding  surface). 


Table  1.  Parameters  used  to  estimate  the  at-sensor  radianee  from  a  seene  eontaining  a  gas 

plume. 


Gas(es) 

Methanol 

Acetone,  Methane,  Sulfur 
Dioxide,  Tetrachloroethylene 

Gas  Emission  Rate(s)  -  g  s'^ 

10 

7,  3,  25,  15 

Gas  Temperature  (C) 

40 

20 

Stack  Height  (m) 

20 

20 

Stack  Diameter  (m) 

0.3 

0.3 

Stack  Exit  Velocity  (m  s'^) 

10 

14 

Atmospheric  Model 

Mid-latitude  Summer 

Sub-Arctic  Summer 

Wind  Speed  (m  s'^) 

5 

10 

Surface  Temperature  (C) 

21.2 

1.2 

Sensor 

Prism  Spectrometer 

FTS 

Sensor  Altitude  (km) 

2.8 

3.5 

Sensor  OPD  (cm) 

- 

0.25 

A  seeond  seenario  eonsisted  of  a  plume  eontaining  four  gases  eaeh  with  a  unique  emission 
rate  as  shown  in  Table  1.  For  this  ease  the  at-sensor  radianee  from  the  plume  and  baekground 
was  modeled  for  a  notional  Fourier  Transform  Speetrometer  (FTS)  with  an  Optieal  Path 
Differenee  (OPD)  of  0.25  em.  This  eorresponds  to  a  speetral  resolution  of  approximately  2 
cm"'.  The  graphs  in  Fig.  4  eorresponding  to  the  derived  TbS,  plume  emissivity  and  absorbance 
are  analogous  to  those  presented  in  Fig.  3.  However,  in  this  case,  overlap  of  the  effluent 
signatures  makes  it  more  difficult  to  identify  individual  contributions.  Fig.  4d  displays  the 
total  plume  absorbance  overlaid  on  top  of  the  individual  absorbances  for  the  four  gases  used 
in  the  simulation.  The  displayed  integrated  concentrations  (in  PPM-m)  were  found  by 
manually  adjusting  the  standard  gas  concentrations  so  that  specific  gas  absorbance  features 
were  in  agreement  with  the  total  plume  absorbance.  Automation  of  this  step  is  under 
development  using  a  spectrally  matched  filter  technique. 

4.  SUMMARY 

Measurements  of  outgoing  radiance  in  the  LWIR  are  known  to  be  useful  for  detection  and 
quantification  of  gaseous  absorbers.  A  model  originally  developed  to  provide  hyperspectral 
sensor  system  analyses  has  been  modified  to  allow  simulation  and  analysis  in  the  LWIR  of 
scenes  containing  gaseous  effluents  emitted  from  a  known  source.  This  model  provides  a 
testbed  to  simulate  the  at-sensor  radiance  measurements  from  current  and  future  sensors  for 
cases  of  plumes  containing  one  or  more  gaseous  absorbers  each  with  unique  emission  rates. 
The  model  computes  the  integrated  effluent  concentrations  and  merges  that  information  with 
atmospheric  optical  depths  computed  from  a  line-by-line  radiative  transfer  model.  A 
component  radiative  transfer  model  allows  estimation  of  the  radiance  components  of  the 
plume,  surface,  atmosphere  and  the  total  at-sensor  amount. 


Equations  describing  the  individual  radiance  components  were  presented  for  transfer  through 
the  plume  and  for  a  clear  atmosphere.  The  differential  radiance  was  computed  and  a 
simplified  form  was  derived  which  allowed  estimation  of  the  contributions  from  the  three 
major  components  to  be  compared;  the  atmospheric  transmittance,  the  plume  emissivity,  and 
the  plume/surface  radiance  difference.  While  the  transmittance  and  plume  emissivity 
contributed  to  the  magnitude  of  the  plume  signal,  the  plume/surface  temperature  difference 
contributed  not  only  to  the  magnitude  but  determined  whether  an  emission  (positive 
differential  radiance)  or  absorption  (negative  differential  radiance)  effect  was  dominant. 

Examples  were  shown  for  cases  with  single  and  multiple  gas  plumes.  Eor  single  gases  of 
sufficient  concentration  and  temperature  difference,  it  is  generally  easy  to  separate  the  plume 
signature  from  the  background  radiance.  Eor  plume  cases  with  multiple  gas  components, 
overlap  of  the  gaseous  absorbance  can  mask  individual  signatures  making  identification  of 


Wavelength  (pm)  Wavelength  (pm) 


Eigure  3.  Plots  of  various  components  as  a  function  of  wavelength  for  a  scene  containing  a 
gaseous  effluent  (Methanol).  The  four  plots  show  a)  the  brightness  temperature  for  both  the 
plume  (solid)  and  background  (dotted),  b)  the  differential  brightness  temperature,  c)  the 
plume  emissivity  and  d)  the  derived  (solid)  and  input  “true”  (dotted)  effluent  absorbance. 


gases  difficult.  A  good  example  of  this  would  be  volcanic  emissions  which  quite  often 
contain  not  only  sulfur  dioxide  (a  primary  constituent  of  volcanic  gas),  but  enhanced  amounts 
of  water  vapor,  carbon  dioxide,  hydrogen  sulfide  and  other  gases  which  are  generally  emitted 
at  temperatures  much  higher  than  the  ambient.  The  broad  extent  of  water  vapor  absorption  in 
the  LWIR  can  make  estimation  of  the  amount  of  the  other  constituents  a  complex 
undertaking.  While  a  manual  technique  was  used  here  to  extract  the  individual  gas  signatures 
from  the  total  plume  absorbance,  other  techniques  such  as  spectral  matched  filter  might 
provide  a  way  to  automate  the  process. 


Figure  4.  Similar  to  Fig.  3  except  the  scenario  includes  four  gases:  acetone,  sulfur  dioxide, 

methane  and  tetrachloroethylene. 
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